HU-EP-13/07, ITEP-LAT/2013-1 



Landau gauge ghost propagator and running coupling in SU{2) lattice gauge 

theory 

V. G. Bornyakov 
Institute for High Energy Physics, 142281, Protvino, Russia 
and Institute of Theoretical and Experimental Physics, 117259 Moscow, Russia 



o 

(N 
X) 

(N 



> 
(N 

o 



E.-M. Ilgenfritz 

Humboldt- Universitdt zu Berlin, Institut fiir Physik, 12489 Berlin, Germany and 
Joint Institute for Nuclear Research, VBLHEP, 141980 Dubna, Russia 

C. Litwinski and M. Miiller-Preussker 
Humboldt-Universitdt zu Berlin, Institut fiir Physik, 12489 Berlin, Germany 

V. K. Mitrjushkin 

Joint Institute for Nuclear Research, BLTP 141980 Dubna, Russia 
and Institute of Theoretical and Experimental Physics, 117259 Moscow, Russia 
(Dated: February 22, 2013) 

We study finite (physical) volume and scaling violation effects of Landau gauge ghost propagator 
as well as of the running coupling as{p) in the SU{2) lattice gauge theory. We consider lattices with 
physical sizes between 3 and 7 fm and values of spacing a between 0.2 and 0.07 fm. To fix the gauge 
we apply an efficient gauge fixing method aimed at finding extrema as close as possible to the global 
maximum of the gauge functional. We show that finite volume effects turn out to be rather small 
for momenta |p| ^ 0.6 GeV. The relative scaling violations for as{p) do not exceed 10% within the 
bare coupling range j3 = G [2.30, 2.55] for the same momentum region and decreases quickly 
with increasing \p\. Our results testify in favor of the decoupling behavior. 
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I. INTRODUCTION 

The infrared (IR) behavior of Landau gauge gluon 
and ghost propagators is believed to be closely related 
to gluon and quark confinement scenarios. For exam- 
ple, the celebrated Gribov-Zwanziger/Kugo-Ojima 
(GZKO) color confinement scenario [ll-Q prescribes 
that the gluon propagator D{p) should vanish in 
the IR limit p (so-called infrared suppression), 
while the ghost dressing function is expected to be- 
come singular (infrared enhancement). In a very re- 
cent paper D. Zwanziger has derived a strict bound 
\imp^op'^~^D{p) = (where d denotes the space-time 
dimension for the Landau gauge case) allowing also 
D{0) ^ for d > 2 [i). 

The search for gluon and ghost propagator solutions 
of Dyson-Schwinger (DS) and functional renormaliza- 
tion group (FRG) equations showed the existence of 
infrared solutions exhibiting a power-like scaling be- 
havior (rl-fisl as well as regular so-called decoupling so- 
lutions providing an IR-finite limit of both the gluon 
propagator and the ghost dressing function |16l - [20l |. 
Both kinds of solutions can be realized by different IR 
boundary conditions for the ghost dressing function as 
it has been argued in [2lj , and both of them can sup- 



port quark confinement [22'|. The decoupling solution 
does not agree with the GZKO scenario. However, 
it is in agreement with the refined Gribov-Zwanziger 
formalism 

From the phenomenological point of view the prop- 
agators can serve as input to bound state equations 
as there are Bethe-Salpeter or Faddeev equations for 
hadron phenomenology d, [1^ [1^ . In the ultraviolet 
limit they allow a determination of phenomenogically 
relevant parameters such as ^jfs '^^ condensates 



e.g. by fitting ab-initio lattice data 



to continuum expressions (see e.g. [27|, |28| and cita- 
tions therein) obtained from operator product expan- 
sion and perturbation theory [29, 30]. 

What concerns the solution of DS and/or FRG 
equations it is well-known that the system of those 
equations has to be truncated. The details of trun- 
cation influence the behavior of the Green functions 
especially in the non-perturbative momentum range 
around 1 GeV, where the Landau gauge gluon dress- 
ing function exhibits a pronounced maximum. There- 
fore, reliable results from ab-initio lattice computa- 
tions to compare with or even used as an input for DS 
or FRG equations are highly welcome. 

On the lattice, over almost twenty years extensive 
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studies of the Landau gauge ghost propagators have 
been carried out (see, e.g., j3ll - l46j ). A serious problem 
in these calculations represent the Gribov copy influ- 
ence [H, El] . Numerical lattice results clearly support 
the decoupling-type of solutions in the IR limit and 
the lack of IR enhancement of the ghost propagator 
p3 - l46l |. However, most of the lattice computations 
dealing with the IR limit were relying on rather coarse 
lattices in order to reach large enough volumes. A sys- 
tematic investigation of lattice discretization artifacts 
or scaling violations was missing for long time. 

In this paper we present an investigation for the 
ghost dressing function and - employing previous 
gluon propagator results (47| - obtain the running 
coupling within the so-called minimal MOM scheme 
(isj . We shall use the same lattice field configura- 
tions as in j47j which were gauge fixed with an im- 
proved method taking into account many copies over 
all Z(2) Polyakov loop sectors and applying simulated 
annealing with subsequent overrelaxation. We sep- 
arately discuss the case of fixed lattice spacing and 
varying volume (from 3 fm to 7 fm) and the case of 
fixed physical volume and varying lattice spacing (be- 
tween 0.21 fm and 0.07 fm). In the range of achieved 
IR momenta finite-size effects are shown to be neg- 
ligibly small. But relative finite-discretization effects 
in the infrared (for a renormalization scale chosen at 
ji = 2.2 GeV) turn out to be more sizable and can 
be quantified to reach a 10 percent variation level at 
p ~ 0.4 GeV in the approximate scaling region ex- 
plored between f3 — 2.3 and /3 ~ 2.55. A similar 
conclusion will be drawn for the running coupling. 

In Section [n] we introduce the lattice Landau gauge 
and the corresponding Faddeev-Popov operator and 
the ghost propagator. In Section IIIII some details of 
the simulation and of the improved gauge fixing are 
repeatedly given. In Section lTVl we present our numer- 
ical results for the ghost propagator and the running 
coupling. Conclusions will be drawn in Section IVl 



II. LATTICE LANDAU GAUGE AND THE 
GHOST PROPAGATOR 

Let us briefly recall how the SU{2) gauge field con- 
figurations used in Ref. [47i| for measuring the gluon 
propagator have been created and gauge fixed. 

The non-gauge-fixed SU{2) gauge field configura- 
tions were generated with a standard Monte Carlo 
routine using the standard plaquette Wilson action 



X fi>f 



(1) 



where go denotes the bare coupling constant. The link 
variables Uxfj. G SU{2) transform under local gauge 



transformations gx as follows 

U., ^ U^^ = glUx,gx+i. , g. e SU{2) . (2) 

The standard (linear) definition [4^ for the dimen- 
sionless lattice gauge vector potential Ax+fi/2.ti is 



A. 



^^{Ux,-Ul)^Al^^,,./-^. (3) 



This definition, which is not unique, can influence the 
propagator results in the IR region, where the contin- 
uum limit is harder to control. 

In lattice gauge theory the most natural choice of 
the Landau gauge condition is by transversality [i^ 

4 

{dA)x = (A+A/2;p - A-A/2;m) = ' (4) 

which is equivalent to finding a local extremum of the 
gauge functional 



(5) 



with respect to gauge transformations gx ■ V = L'^ 
denotes the Ad lattice size. The Gribov ambiguity 
is reflected by the existence of multiple local extrema. 
The manifold consisting of Gribov copies providing lo- 
cal maxima of the functional ([5]) and a semi-positive 
Faddeev-Popov operator (see below) is called the Gri- 
bov region Q, while the global maxima form what is 
called the fundamental modular region (FMR) A C 
Our gauge fixing procedure is aimed to approach A by 
finding higher and higher maxima. This is achieved by 
use of the effective optimization algorithm and finding 
a large number of local maxima of which the highest 
is picked up. For recent alternative approaches see 
|50,IIll- 

The lattice expression of the Faddeev-Popov opera- 
tor M"'' corresponding to M"'' = —9^1?^'' in the con- 
tinuum theory (where D'^'' is the covariant derivative 
in the adjoint representation) is given by 



j^rab _ r / qab , cab 

xy / J [ V a^M x — (i; 



x+fi 



(6) 



where 



ST, = 5"" - Tr Ux,, 



Aah ^ 



abc 



a; 



x-'rfi/2:fi' 



(7) 



From the form ([7]) it follows that a trivial zero eigen- 
value is always present, such that at the Gribov hori- 
zon dT the first non-trivial zero eigenvalue appears. 
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For configurations with a constant field, witli 6°^ = 
and 6°^ = 6° independent of x, tliere exist eigenmodes 
of M witli a vanislring eigenvalue. Thus, if the Landau 
gauge is properly implemented, -M[C/] is a symmetric 
and semi-positive definite matrix. 

The ghost propagator G°-^{x, y) is defined as [biIIs^I 



G'^\x,y) = S^' G{x - y) ^ ( (M^^ )l\ [U] 



(8) 



where M\U] is the Faddeev- Popov operator, on the 
sector orthogonal to the strict zero modes. Note that 
the ghost propagator becomes translationally invari- 
ant (i.e., dependent only on x — y) and diagonal in 
color space only in the result of averaging over the 
ensemble of gauge-fixed representants of the origi- 
nal gauge-unfixed Monte Carlo gauge ensemble. The 
ghost propagator in momentum space can be written 
as 



X, V 



"'"^^■^"-^)((Af-i)^;[C/]) , (9) 



where the coefficient is taken for a full normal- 
ization, including the indicated color average over 
a = 1, ..,3. In what follows we will denote the (bare) 
ghost dressing function as 



J{p)=p'Gip). 



(10) 



III. DETAILS OF THE COMPUTATION 

The Monte Carlo (MC) simulations had been car- 
ried out at several /3-values between /3 = 2.2 and 
(3 = 2.55 for various lattice sizes L. Consecutive 
configurations (considered to be statistically indepen- 
dent) were separated by 100 sweeps, each sweep con- 
sisting of one local heatbath update followed by L/2 
microcanonical updates. In TableUwe provide the full 
information about the field ensembles used in this in- 
vestigation. The corresponding information concern- 
ingj ^ the gluon propagator has been communicated in 

we employ the Z{2) fiip operation 
_ [5^. For completeness we recall 
The method consists in flipping all 
attached and orthogonal to a se- 
lected 3d plane by multiplying them with —1 £ Z(2). 
Such global flips are equivalent to non-periodic gauge 
transformations. They represent an exact symmetry 
of the pure gauge action. The Polyakov loops in the 
direction of the chosen links and averaged over the or- 
thogonal 3d plane obviously change their sign. There- 
fore, the fiip operations combine the 2^ distinct gauge 
orbits (or Polyakov loop sectors) related to strictly pe- 
riodic gauge transformations into a single large gauge 
orbit. 



For gauge fixing 
as discussed in [5£ 
the main features, 
link variables Uxf_i 



p 


[GeV] 


a [fm] 


L 


aL [fm] 


^meas 


copy 


2.20 


0.938 


0.210 


14 


2.94 


400 


48 


2.30 


1.192 


0.165 


18 


2.97 


200 


48 


2.40 


1.654 


0.119 


26 


3.09 


200 


48 


2.50 


2.310 


0.085 


36 


3.06 


400 


80 


2.55 


2.767 


0.071 


42 


2.98 


200 


80 


2.20 


0.938 


0.210 


24 


5.04 


400 


48 


2.30 


1.192 


0.165 


30 


4.95 


400 


48 


2.40 


1.654 


0.119 


42 


5.00 


200 


80 


2.30 


1.192 


0.165 


44 


7.26 


200 


80 



TABLE I: Values of /3, lattice sizes, number of measure- 
ments and number of gauge copies used throughout Ref. 
[47I ] and this paper. The lattice spacing was fixed to its 
physical value using the string tension ^/a = 440 MeV (see 
[53,,^]). 



The second ingredient is the simulated annealing 
(SA) method, which has been investigated indepen- 
dently and found computationally more efficient than 
the exclusive use of standard overrelaxation (OR) (sgI - 
[ssl ]. The SA algorithm generates gauge transforma- 
tions g{x) by MC iterations with a statistical weight 
proportional to eyi\){W Fu[g\/T) . The "tempera- 
ture" T is an auxiliary parameter which is gradu- 
ally decreased in order to guide the gauge functional 
Fjj[g\ towards a maximum, despite its fiuctuations. 
In the beginning, T has to be chosen sufficiently large 
in order to allow rapidly traversing the configuration 
space of g{x) fields in large steps. As in Ref. [56] we 
have chosen Tinn = 1.5. After each quasi-equilibrium 
sweep (that includes both heatbath and microcanon- 
ical updates) T has been decreased in equidistant 
steps. The final SA temperature has been fixed ac- 
cording to the requirement that during the following 
execution of the OR algorithm the violation of the 
transversality condition 



max 

a:, a 



4 

E 



+A/2;m 



-/i/2;/x 



(11) 



decreases in a monotonous manner for the majority 
of gauge fixing trials, until finally the transversality 
condition (|lip becomes uniformly satisfied with an 
eior = 10~^. A monotonous OR behavior is reason- 
ably satisfied for a final lower SA temperature value 
T final = 0.01 57]. The number of temperature steps 
interpolating between Tina and T final has been cho- 
sen to be 1000 for the smaller lattice sizes and in- 
creased to 2000 for the lattice sizes 30^ and bigger. 
The finalizing OR algorithm using the Los Alamos 
type overrelaxation with the overrelaxation parameter 
value a; = 1.7 requires typically a number of iterations 
varying from 0(10^) to 0(10^). 
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In what follows we call the combined algorithm em- 
ploying SA (with finalizing OR) and Z{2) flips the 
'FSA' algorithm. By repeated starts of the FSA al- 
gorithm we explore each Z{2) Polyakov loop sector 
several times in order to find there the best ("6c ") 
copy ^. The total number of copies per configuration 
Ncopy, generated and inspected for each /3- value and 
lattice size is indicated in Table HI 

Some more details suitable to speed up the gauge 
fixing procedure are described in j45j. 

In order to suppress lattice artifacts from the be- 
ginning we followed Ref. [s^ and selected the allowed 
lattice momenta as surviving the cylinder cut 



the data &t (3 — 2.3 for aL ~ 7 fm and at /3 = 2.4 for 
aL ~ 5 fm with a fitting function of the form 



(12) 



Moreover, we have applied the "a-cut" [59| < 
(2/a)a for every component, in order to keep close 
to a linear behavior of the lattice momenta = 
(27rfc^)/(aL), k^, e (-i/2,L/2]. We have chosen 
a = 0.5. Obviously, this cut infiuences large momenta 
only. 

We define the renormalized ghost dressing function 
according to momentum subtraction schemes (MOM) 
by 

Jrenip, m) = Z {fl, 1 / a) J {p, I / a) (13) 

Jrenip = ^l) = 1. (14) 

In practice, we have fitted the bare dressing function 
J{p,\/a) with an appropriate function (see Eq. (|15p 
below) and then used the fits for renormalizing J. 
Assuming that lattice artifacts are sufficiently sup- 
pressed it has to be seen, whether multiplicative 
renormalizability really holds in the non-perturbative 
regime. For this it is sufficient to prove that ratios 
of the renormalized (or unrenormalized) propagators 
obtained from different cutoff values l/a(/3) will not 
depend on p at least within a certain momentum inter- 
val [pmin-,Pmax\-, whcrc Pmax should bc the maximal 
momentum surviving all the cuts applied. 

In what follows the subtraction momentum has 
been chosen to be ji — 2.2 GeV. 



IV. RESULTS 
A. Finite volume effects 

Data for Jren (p, A*) for various volumes are pre- 
sented in Fig. [T] for P = 2.3 and in Fig. H for /3 = 2.4. 
To see the finite volume effects in more detail we fitted 



fj{p) 



hi 



(p2 



(15) 



gh 



^ An alternative idea to tackle the Gribov problem has been 
discussed in ISOl. ISlll. 



(see Table HI)) . This ansatz, while describing the data 
reasonably well within the given momentum range, 
will not be applicable in the IR limit, when we as- 
sume that J(p) exhibits an inflection point and bends 
to a finite value J(0). The latter was reported for 
SU{2) [45] at ^ = 2.20 on larger lattice sizes (40'') 
than considered here. 

In the right panels of Fig. [1] and Fig. [21 respectively, 
the relative deviations from the fit function are shown 
for (i = 2.3 and /? = 2.4, respectively. One can see 
that for both /3 values finite volume effects for lattices 
even with ai ~ 3 fm are small (less than 1%) for 
momenta \p\ ^ 0.6 GeV. 



B. Lattice spacing effects 

In Fig.[3]we show the momentum dependence of the 
renormalized ghost dressing function Jrenip) for five 
different lattice spacings but for (approximately) the 
same physical size aL ~ 3 fm (for the exact values see 
Table U. Finite-spacing effects for /? = 2.2,2.3,2.4 in 
comparison with /3 — 2.55 are evident. To see them in 
more detail we have fitted also the data at /3 = 2.55 
with the fit function Eq. (see Table HH). Then, in 
the right panel of Fig. |3] we have plotted the relative 
deviation of the data from this fitting result for all 
values of the lattice coupling parameter /3. Note that 
the deviations of the data points for (3 = 2.55 from 
the zero-constant line just illustrate the quality of the 
fit. 

Related to our choice of the (re)normalization mo- 
mentum ^ = 2.2 GeV and due to the rather small 
statistical errors for the ghost dressing function we 
see clear scaling violations especially in the IR region 
but also for p > /i. At the lowest (here accessible) 
momenta the violations a.t /3 — 2.3 {(3 = 2.4) rela- 
tive to ^ = 2.55 are staying below 10% (4%). Thus, 
in comparison with corresponding estimates for the 
gluon propagator (see Fig. 13 in which were 

more noisy, we can say that the relative scaling vi- 
olations of the ghost dressing function turn out to be 
somewhat larger. 

Similar to the case aL ~ 3 fm we observe analogous 
lattice spacing effects on volumes with linear size aL ~ 
5 fm. The respective results are depicted in Fig. |4l 
In this case the fit for (3 — 2.4 was employed (see 
Tabic HH. 
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FIG. 1: Left: The momentum dependence of the renormalized ghost dressing function Jren{p) for three different lattice 
sizes at P = 2.3. The curve shows the fit applying Eq. (|f 5p to the case aL ~ 7 fm. Right: The relative deviation of the 
data for the ghost dressing function Jren(p) from the former fitting curve. 




/? 


L 


■nigh 




bi 


b2 


Xdf 


2.55 


42 


0.67(1) 


0.0211(2) 


1.92(2) 


-0.94(2) 


0.99 


2.40 


42 


0.67(1) 


0.0232(4) 


2.01(1) 


-1.03(1) 


1.9 


2.30 


44 


0.64(1) 


0.026(1) 


2.15(1) 


-1.15(1) 


1.8 



TABLE IL Values of the fit parameters in Eq. (fT5|l in 
physical units and the corresponding Xdf- 



C. Running coupling 

Taking the gluon dressing function results from 
[47I I into account we can compute the minimal MOM 
scheme running coupling psl j via 



(16) 



where Z(p) and J{p) are the bare gluon and ghost 
dressing functions, respectively. 

For the running coupling we use the following fitting 
function: 



Up) 



cip 



C2P 



C3P 



ml: 



{p'^+m^y ' (p2-|-m^)3 



(17) 



The fit results for the same combinations of values 
{(3, L) as for the ghost dressing function (see Table HI)) 
are collected in Table Hill 

Finite volume effects for the running coupling are 
shown in Fig. [S] for (3 — 2.3 and in Fig. [5] for (3 = 2.4. 
In both cases one can see the finite volume effects to 
be reasonably small (less than 5%) at a linear physical 
lattice extension aL ~ 3 fm and for momenta \p\ ^ 0.6 
GeV. 

Results for the scaling check of as{p) taking into 
account four lattice spacings for the same extent of 
aL ~ 3 fm are presented in Fig. [71 We see relative 
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FIG. 3: Left: The momentum dependence of the renormalized ghost dressing function Jrenijp) for five different /?- values 
or lattice spacings. The physical linear box size is aL ~ 3 fm. The fitting curve belongs to the smallest available lattice 
spacing (/? = 2.55). Right: The relative deviation of the ghost dressing function data Jren(p) from the fitting result. 




0.5 1 1.5 2 2.5 3 3.5 0.5 1 1.5 2 2.5 3 3.5 

p, GeV p, GeV 

FIG. 4: Same as in Fig. |3]but for a linear box size of aL ~ 5 fm and three different /3- values. The fitting curve belongs 
to the smallest available lattice spacing (/3 = 2.40). 



deviations for /3 = 2.3 in comparison with /3 = 2.55 up 
to a 10%-level within the momentum range explored. 

Since the running coupling as seems to tend to 
zero in the IR limit, our results obtained within the 
framework of Landau gauge fixing as described above 
are fully compatible with the IR-decoupling scenario 
discussed in the context of the Dyson-Schwinger and 
functional renormalization group approach (2ll [60j . 





L 


ma 


Cl 


C2 


C3 


2 
Xdf 


2.55 


42 


1.04(1) 


0.205(3) 


2.47(3) 


13.9(9) 


1.3 


2.40 


42 


1.01(1) 


0.16(1) 


2.68(5) 


11.7(7) 


0.81 


2.30 


44 


1.03(1) 


0.19(4) 


2.4(2) 


14(2) 


0.95 



TABLE III: Values of the fit parameters for Us (Eq. (fTT)) ') 
in physical units and the corresponding Xd/ • 



V. CONCLUSIONS 

Completing an earlier work [l^] we have computed 
the Landau gauge ghost dressing function for lattice 
SU (2) pure gauge theory. Together with the former 
results for the gluon propagator we also have pre- 
sented the minimal MOM scheme running coupling. 
We employed the same gauge-fixed field configura- 
tions as in [i^] obtained with a gauge fixing method 
consisting of a combined application of Z(2)-flips and 
repeated simulated annealing with subsequent overre- 
laxation for the gauge functional. This method was 
invented in order to get as close as possible to the fun- 
damental modular region i.e. to the global extremum 
of the gauge functional. It previously was shown to 
solve the Gribov problem with suppressed finite size 
effects [56l, Issj. 

Assuming that the Gribov problem is kept under 
control to the best of our present knowledge we con- 
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0.5 1 1.5 2 2.5 0.5 1 1.5 2 2.5 

p, GeV p, GeV 

FIG. 5: Left: The momentum dependence of the running coupling for three different lattice sizes at /3 = 2.3. The curve 
shows a fit with Eq. (jf 7[l for aL ~ 7 fm. Right: The relative deviation of the data for the running coupling aa{p) from 
the fitting curve. 




centrated ourselves on systematic effects like finite size 
and lattice spacing dependences. While finite size ef- 
fects were confirmed to be rather small the lattice 
spacing artifacts turned out to be non-negligible as 
for the renormalized ghost dressing function as for 
the running coupling. In both cases for the linear 
lattice size of approximately 3 fm (and for the ghost 
dressing function with a subtraction momentum of 
/i = 2.2 GeV) we saw relative deviations at /3 = 2.30 
from the results obtained at our largest /3 = 2.55 
reaching a ten percent level at the lowest accessible 
momentum values and still around five percent in the 
non-perturbative region around 1 GeV. This tells us 
that lattice results for Landau gauge gluon and ghost 
propagators in this momentum range have still to be 
taken with some caution what concerns the Gribov 
problem and the continuum limit. 

However, even if the ghost dressing function for not 
too small momenta has been well fitted by a weakly 



IR singular behavior (see Eq. (|15p ). thus hiding a IR 
finite limit, the result for the running couling as{p) 
turned out to be robust what concerns the compati- 
bility with the infrared decoupling solution of DS and 
FRG equations. 
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012345 012345 
p, GeV p, GeV 

FIG. 7: Left: The momentum dependence of the running couphng for four different /3- values or lattice spacings at lattice 
extent aL ~ 3 fm. The curve shows the fit result with Eq. (jlTp for /3 = 2.55. Right: The relative deviation of the data 
from this fit. 
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